Aging and multiscaling in out of equilibrium dynamical processes of granular media 
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In the framework of recently introduced frustrated lattice gas models, we study the out of equilib- 
rium dynamical processes during the compaction process in granular media. We find irreversible- 
reversible cycles in agreement with recent experimental observations. Moreover in analogy with 
the phenomenology of the glass transition we find aging effects during the compaction process In 
particular we hnd that the two time density correlation function C(t,t') asymptotically scales as 
a function of the single variable ln(t')/ ln(t). This result is interpreted in terms of multiscaling 
properties of the system. 



The experimental study of dynamic processes in gran- 
ular media Q has recently revealed the presence of in- 
teresting behaviours. Under tapping dry granular media 
reach very slowly a more compact state which is well fit- 
ted by a logarithmic relaxation Q . More recently Novak 
et al. H have also shown that such materials exhibit non 
trivial irreversible-reversible cycles. These phenomena 
stem from slow relaxation processes due to large "co- 
operative rearrangements" of many particles. In such a 
perspective granular materials share features of thermal 
systems such as glasses or spin glasses which are also 
characterized by extremely long relaxation times and the 
presence of irreversible-reversible cycles . 

In this paper in the framework of recently introduced 
microscopic models we reproduce the irreversible- 

reversible cycles of Novak et al. and investigate the effect 
of the "cooling" rate on the compaction process. We find 
a behaviour which is strongly reminiscent of the phe- 
nomenology of the glass transition. Finally we study the 
non equilibrium time dependent density density autocor- 
relation function and found aging effects typical of glassy 
systems. These results suggest that similar effects could 
be also found in real experiments. 

In dynamical processes of granular media a crucial role 
is played by geometric frustration (originated by steric 
hindrance between interlocked neighboring grains) which 
induce the necessity of large scale cooperative rearrange- 
ments for relaxation. Based on these concepts two kinds 
of frustrated lattice gas model were introduced ji] |). 
Both models showed logarithmic compaction, segregation 
and other phenomena typical of granular media under 
shaking. 

These models consist of a system of particles which 
occupy the sites of a square lattice tilted by 45° (see 
Fig.|l|). Particles are characterized by an internal degree 
of freedom, Si = ±1, corresponding for instance to two 
typical orientations of grains on the lattice. Two nearest 
neighbor sites can be both occupied only if the parti- 
cles have the right reciprocal orientation, so that they 
do not overlap, otherwise, due to excluded volume, they 
have to move away. In absence of vibrations the particles 
are subject only to gravity and they move downwards al- 
ways fulfilling the non overlap constraint. The effect of 



vibration is introduced by allowing the particles to dif- 
fuse with a probability p up upwards and a probability 
Pdown = 1 — Pup downwards. An important parameter 
governing the dynamics is the adimensional parameter 
r = — l/ln(a;o), with xq = p U p/pdown, which is related 
(see below) to the effective temperature of the system 
and consequently plays the same role as the amplitude of 
the vibration intensity in the experiment of Novak et al. 




FIG. I. A schematic picture of the two kind of frus- 
trated lattice gas models described in the text. Left: the 
Tetris model. Right: the Ising frustrated lattice gas, IFLG. 
Straight and dashed lines represent the two kind of inter- 
actions €ij — ±1. Filled circles are present particles with 
"orientation" Si — ±1 (black/white). 

Such models can be described in terms of the following 
lattice gas Hamiltonian (see ^,|)) in the limit J — > oo: 

H = J );fij{Si,Sj)niiij (I) 

(ij) 

Here ni = 0, f are occupancy variables, Si = ±1 are 
spin variables, associated to the two orientations of the 
particles, J represents the infinite repulsion felt by the 
particles when they have the wrong orientations and 
fij(Si, Sj) = or 1 depending whether the configuration 
Si,Sj is right (allowed) or wrong (not allowed). 

The choice of fiASi,Sj) depends on the particular 
model. The Tetris g model is made of elongated parti- 
cles (see Fig. Q), which may point in two directions coin- 
ciding with the two lattice bond orientations. In this case 
fij(Si,Sj) is given by fij(Si,Sj) = l/2(S l S J - e i . j (S l + 
Sj) + 1), where = +1 for bonds along one direction 
of the lattice and e,j = — 1 for bonds on the other. This 
hamiltonian model has an ordered "antiferromagnctic" 
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ground state, and its dynamics has the crucial constraint 
that particles can flip their "spin" only if 3 of their own 
neighbors are empty. 

A real granular system may contain more disorder due 
to a wider shape distribution or to the absence of a lat- 
tice. Therefore in a more realistic model the number of 
internal states is q > 2 (Si = l,2...q) and the function 
fij(Si, Sj) is zero only for allowed nearest neighbor con- 
figurations. 

However to simplify the model the number of states 
was still kept q = 2 and the randomness was taken into 
account by introducing random quenched variables corre- 
sponding to the freezing of some degree of freedom in the 
high density regime. Thus an Ising frustrated lattice gas 
model (IFLG) was proposed Jt| in which fij(Si, Sj) was 
given by fij(Si,Sj) = l/2(e lj S l S 3 - 1) and ey = ±1 are 
quenched random interactions associated to the bonds of 
the lattice. 

The Hamiltonian (0) is without gravity. In presence 
of gravity there is an extra term g ^ riiUi where g is 
the gravity and yi is the ordinate of the particle i. The 
temperature T is related to the ratio xq = p U p/Pdown via 

e -2g/T = ( notice p = T / 2 g ). 

Interesting enough the two extreme models, the Tetris 
and the IFLG, show similar behaviour. This suggests 
that the results found are rather robust and will not de- 
pend much on the details of the model. In particular 
under tapping they reproduce the logarithmic behaviour 
in agreement with the experimental results of Knight et 
al. j^] . Experimentally a "tap" is the shaking of the con- 
tainer of the grains by vibrations of given duration and 
amplitude. In our Monte Carlo simulations each single 
tap is realized by letting the particle diffuse under the 
gravity by keeping T = const, during the time interval of 
a tap, To, and then switching off the vibration by setting 
r = until the system reaches a static configuration. 
Time t is measured in such a way that one unit corre- 
sponds to one single average update of all particles and 
all spins of the lattice. 

In a previous paper Q we have performed on the IFLG 
a particular cycle sequence of taps to show the presence in 
granular media of hysteresis effect. Recently also real ex- 
periments were performed on density relaxation in gran- 
ular media under cyclic tapping ||], which show indeed 
the presence of such hysteresis effect. In order to com- 
pare better with the experimental data we have simulated 
the same cycle of Novak et al. on the IFLG model. We 
have considered a system of size 30 x 60 (our data are 
robust to size changes), with periodic boundary condi- 
tions along the x-axis and rigid walls at bottom and top. 
Our data are averaged over 8 different lattice realizations 
each averaged over 30 different noise realizations. 

A starting particle configuration is prepared by ran- 
domly inserting particles into the box from its top and 
then letting them fall down, with the dynamics described 
above, until the box is filled. We performed cycles of taps 
in which the vibration amplitude T is varied at fixed am- 



plitude increment 7 = AT/tq holding constant their du- 
ration To. More precisely we performed a sequence of J\f 
taps, of amplitude ri...F„...rjv, from an initial ampli- 
tude Ti = to a maximal amplitude T max = 15, then 
back to r = and then up again to Tjv = T max . After 
each tap we measure the static bulk density of the system 
p(T n ) (n is the n-th tap number). 
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FIG. 2. The static bulk density, p(T), of the IFLG model 
(the Tetris gives analogous results) as a function of the vi- 
bration amplitude, T, in cyclic vibration sequences. The 
system is shaken with an amplitude V which at first is in- 
creased (filled circles), then is decreased (empty circles) and, 
finally, increased again (filled squares) with a given "anneal- 
ing-cooling" velocity 7 = AT /to (at each value of F the sys- 
tem has undergone a "tap" of duration To = 10 3 ). Here we 
fixed 7 = 1.25 10 -3 . The upper part of the cycle is approxi- 
mately "reversible" (i.e., empty circles and filled squares fall 
roughly on the same curve). The data compare rather well 
with the experimental data of Novak et al.. F* is approxi- 
mately the point where the "irreversible" and the "reversible" 
branches meet. r 9 signals the location of a "glass transition" 
(see Fig.|). 

Our results are qualitatively very similar to those re- 
ported in real experiments on dry granular packs j|] . We 
find that when the system is successively shaken at in- 
creasing vibration amplitudes, the bulk density of the 
system typically grows and then decreases as shown in 
Fig. ||. However, when the amplitude of shaking de- 
creases back, the density follows the same path only up 
to some value of T and then deviates from it, in fact it 
does not bend down and keeps growing. As in the ex- 
perimental data the second part of the shaking cycle is 
approximately reversible (see Fig. |^). For such a reason 
these processes are called "irreversible-reversible" cycles 
H . Interestingly the reversible cycle is a monotonic func- 
tion of the shaking amplitude (see FigJ^). 
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To study the dependence on the "cooling" rate 7, we 
repeated the tapping sequence for different values of the 
tap amplitude increment Ar with a fixed tap duration To 
(see Fig.||). We find that the reversible branches have a 
common part for high values of T while for small T they 
splits in different curves depending on the cooling rate 7. 
The slower the cooling rate the higher is the the final den- 
sity observed at the end of the descending part of the cy- 
cle. One can schematically define the point, r g (j), where 
the system freezes and goes out of equilibrium as the lo- 
cation of the "shoulder" in these "reversible" branches 
(see Fig.||and|3|). Thus, 1^(7), which depends on 7, cor- 
responds to a "glass transition". Notice that 1^(7) is 
usually different from the point, T*(7), where the "irre- 
versible" and the "reversible" branches meet (see Fig.|]). 
However, as 7 gets smaller T g and F* become closer and 
they may coincide in the limit 7^0. 
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FIG. 3. As in Fig.H, we report the density, p(T), as a func- 
tion of the vibration amplitude, F for three different values 
of the "cooling" velocity 7. For sake of clarity, we plot here 
only the descending reversible parts of the cycle. As in ex- 
periments on glasses a too fast cooling drives the system out 
of equilibrium. The position of the shoulder, r s (7), in these 
curves schematically individuates a "glass transition" . 

As in glasses the system gets out of equilibrium due 
to the fact that the characteristic times of relaxation are 
much larger than the time To involved in the experiment, 
but, for the same reason, the location of the path depends 
on the rate 7. 

The limit of 7 going to zero defines an ideal glass transi- 
tion amplitude Tq. Using the analogy with the glass tran- 
sition we expect a slow logarithm dependence of T 9 (7) on 
7 as in the glass transition. Full details will be presented 
elsewhere [[LOl . Experimental results in this direction will 
be also very interesting. 

By further exploiting the analogy with the glass tran- 



sition, we expect aging phenomena. For instance if one 
keep shaking the system at a low fixed amplitude T on 
the reversible branch, the system will slowly (logarith- 
mically) approach an asymptotic value of the density, 
eventually on the equilibrium curve. In order to further 
quantitatively characterize the out of equilibrium dynam- 
ics and aging phenomena in granular matter and make 
quantitative predictions, in analogy with glassy system 
we introduce a two time density-density correlation func- 
tion (t > t')\ 



C(t,t>) = 



(p(t)p(t'))-(p(t))(p(t')) 

W?) - W)) 2 



(2) 



where pit) is the bulk density of the system at time t. In 
out of equilibrium C(t,t') is a function of both times t 
and t' (at equilibrium just of t — t'). The aging properties 
of the system are characterized by the specific scaling 
properties of C(t, t 1 ) 

In order to study the system in a well defined config- 
uration of its parameters, we evaluate C(t,t') during a 
"single tap" : we prepare the system at t = by ran- 
domly pouring grains in the box from above as described 
before, then we start to shake it continuously and indefi- 
nitely with a given (small) amplitude T. We expect very 
similar results by considering, instead of a long tap, a 
series of short taps which is experimentally more conve- 
nient (as in ref. g). The data about C(t,t') we present 
here are averaged, at least, over 8 different lattice and 
512 different noise realizations. 

At low r, a good fit for the two time correlation func- 
tion, C(t, t'), on the whole five decades in time explored, 
is given by the following: 



C(M') = (1 - Coo) 



\n{(t> + t s )/r] 
ln[(t + t a )/r] 



(3) 



where r, t s and Coo are fit parameters. Very interesting is 
the fact that the above behavior is found in both of our 
models (Tetris and IFLG). The data for the two models, 
for several values of T, rescaled on a single universal mas- 
ter function, are plotted in Fig. ||. In particular, in the 
explored range of T £ [0.11,0.43] (i.e., x £ [10~ 4 , 10" 1 ]), 
we found r ~ e z / r (with z ~ 2 in both the IFLG and the 
Tetris), the mark of activated dynamics (t s (r) behaves 
approximately as r as a function of T). The asymp- 
totic value Coo is difficult to determine with some pre- 
cision: in the above T range we evaluated approximately 
Coo = 0.2 -r 0.3 for the IFLG model and Coo = 0.0 0.2 
for the Tetris model. 

Eq. (H) essentially states that for times long enough 
the correlation C(i, t') is a function (linear) of the ra- 
tio ln(i')/ ln(i). Such a scaling behaviour is known in 
others disordered systems like Random Ferromagnets or 
Random Fields models |llj] and has been proposed by 
Fisher and Huse droplets theory of finite dimensional 
spin glasses (JT2). However, it seems to be different from 
other scaling functions proposed to fit experimental data 
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in spin glasses M . All this shows the necessity of exper- 
imental confirmation of our results in the framework of 
granular media. 

In any case, eq. (0) is in good agreement with a general 
multiscaling approach presented in Refs. When- 
ever the characteristic relaxation times of the system be- 
comes exceedingly large respect to those typical of the 
measurements, it is reasonable to assume scaling prop- 
erties for the system. The most general scaling trans- 
formation which satisfies group properties gives rise to 
a multiscaling form of the correlation function given by 

& 

C{t,t')=C(a)t f{a) (4) 

with a twofold possibility for a: a = t'/t z or a = 
ln(t')/Ln(£) + 0(l/ln(t)); here C(a) and f(a) are two 
generic functions. The assumption of multiscaling might 
thus naturally allow to explain the presence of large 
classes of universality found in scaling behaviour of ap- 
parently different systems. The above case of eq. (||), 
observed in both Tetris and IFLG, corresponds to a = 
ln(i')/ln(t) and f(a) ~ 0. 
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FIG. 4. The two time density-density correlation function, 
(C(t,t') — Coo)/(l — Coo), as a function of the scaling variable 
a = ln[(t + ts)/r]/ ln[(i' + t s )/r]. Scaled on the same master 
function are data from both models considered in the present 
paper (Tetris, squares, and IFLG, circles) and for vibration 
amplitudes T = — l/ln(xo) with xo € [10 -4 , 10 -1 ]. The mas- 
ter function is 1/q. Inset: The correlation C(t,t') for the 
Tetris at F = 0.22 (or xo ~ 0.01) as a function of t — t' for 
four values of t' = 10 2 , 2 10 2 , 10 3 , 10 4 . 



mic compaction or segregation. Here, we have studied 
irreversible-reversible cycles and found good agreement 
with the experimental data on granular packs ||. We 
have stressed the strong analogies with non equilibrium 
dynamic properties observed in glassy systems. In partic- 
ular we evaluated the two time density correlation func- 
tion C{t,t'), which are found asymptotically to be func- 
tion of the single ratio a = ln(i')/ ln(i). This results ob- 
served in different models, and amenable to experimen- 
tal check, may be interpreted by assuming multiscaling 
properties of the system. 

We thank INFM-CINECA for computer time on Cray- 
T3D/E. 
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In conclusion, in the framework of simple frustrated 
lattice gas models, we have studied the off equilib- 
rium dynamics of slightly shaken granular materials. 
These models have previously shown to share many 
phenomena characteristic of granular media as logarith- 
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